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Abstract 

We consider dilepton production in high energy proton-nucleus (and 
very forward nucleus-nucleus) collisions. Treating the target nucleus as a 
Color Glass Condensate and describing the projectile proton (nucleus) as 
a collection of quarks and gluons as in the parton model, we calculate the 
differential cross section for dilepton production in quark-nucleus scatter- 
ing and show that it is very sensitive to the saturation scale characterizing 
the target nucleus. 



1 Introduction 

Perturbative QCD predicts a sharp rise in the number of gluons per unit area 
and rapidity in a proton/nucleus at high energy (small x) jjj. This however 
would lead to violation of unitarity of hadronic cross sections at high energies. 
High gluon density and gluon recombination effects are believed to be respon- 
sible for taming this growth and restoration of unitarity. It has been suggested 
that gluons at small x can be described by a strong classical field || and that 
weak coupling, semi-classical methods can be applied to describe the physics of 
dense gluonic systems, such as a proton or nucleus at high energies. This dense 
system of gluons (the Color Glass Condensate) is characterized by a saturation 
momentum Q s (x) which grows fast with energy and rapidity. 

There has been much work done in order to investigate the properties of the 
Color Glass Condensate in DIS electron-proton and electron- nucleus, as well as 
proton-nucleus and nucleus- nucleus collisions 0). While the saturation effects 
in protons at current energies are far from being established , the situation in 
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high energy nucleus- nucleus collisions is more intriguing. The saturation model 
seems to work reasonably well at RHIC || even though there are a lot of open 
questions which need to be addressed before one can claim that the saturation 
model describes high energy heavy ion collisions at RHIC quantitatively. The 
role of final state interactions, thermalization, etc. is still to be understood |Q. 

In a recent set of publications j|, ||] , we proposed that high energy proton- 
nucleus collisions at forward rapidities at RHIC may be an ideal place in order 
to investigate the Color Glass Condensate and the saturation model. By con- 
sidering proton-nucleus collisions in the forward rapidity region, one can avoid 
most, if not all, of the complications present in a nucleus- nucleus collision. Pre- 
viously, we calculated (real) photon production rate in p — A collisions || and 
showed that it is very sensitive to saturation effects in the nucleus. 

In this work, we consider dilepton (virtual photon) production in p — A and 
show that dileptons provide a more versatile probe of the saturation model than 
photons. Furthermore, photons are notoriously difficult to measure in a collider 
environment. One needs to define isolation criteria in order to separate photons 
from different sources which greatly reduces the production rates in addition to 
introducing theoretical ambiguities in defining the isolation criteria. The use of 
factorization theorems, well established for highp^ inclusive photon production, 
may be questionable for isolated photons. By considering dilepton production, 
one can avoid most of these experimental and theoretical difficulties [II . 

We briefly review our formalism and the differences between real and virtual 
photon cross sections in Section 2. In section 3, we consider the diffractive 
cross section and show that it vanishes. We consider the inclusive cross section 
in Section 4 and derive the differential cross section for dilepton production in 
p — A collisions da / dz dM 2 d 2 k±, where M 2 and k± are the dilepton invariant 
mass and transverse momentum while z is its fractional energy. We end by 
discussing our results and the experimental signatures of saturation effects. 

2 Real vs. virtual photon production 

In order to reuse some parts of the calculation we already performed for real 
photon production |J, we start by a section highlighting the main differences 
between real photon and lepton pair production, as well as the common aspects. 
We want to calculate the amplitude for the elementary process 

q( P ) + A^q(q)+l+(k 1 )+l-(k 3 )+X (1) 

where a quark entering in the color field of a nucleus emits a lepton pair l + l~ . 
In terms of in- and out- states, this amplitude reads: 

(q(q)l + {k^r (k 2 ) out \q(p) in ) = (Oont\bout(q)bl n (p)c ou t(k 2 )d out (k 1 )\0 in ) , (2) 

where w is the creation operator for a quark, while c 1 and d^ respectively cre- 
ate a lepton and an anti-lepton. Applying the LSZ reduction formula to this 
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Figure 1: A typical contribution to the cross-section for lepton pair production 
in pA collisions. The black dot denotes the re-summed interactions of the in- 
coming quark with the classical color field of the nucleus. This is the square of 
the term where the photon is emitted after the scattering on the nucleus. There 
is also a term where the photon is emitted first, and interferences thereof, not 
represented here. 



amplitude |12[| , one obtains the following expression in terms of the fermionic 
fields: 

(0 O ut|&out(g)^ n (P) c out(fc 2 )rfout(fci)|0in) = J d A xd i yd i z l d A z 2 

xe l ^ x - p - y) e^ kl - Zl+k2 - Z2) u(q)w{k 2 ){i X -m)(i Z2 -m) 

x(0 out |T^(x)^(y)*(z 1 )*(z 2 )|0 in )( l y +m)(i Zl +m)v{k 1 )u{p) , 

(3) 

where tp is the quark field, ^ the leptonic field, u a quark free spinor, w a lepton 
free spinor and v an anti- lepton free spinor. Note that we have approximated 
all the renormalization constants by 1 since we are going to compute only the 
lowest order in the couplings a cm and a s . Therefor this amplitude is made of a 
quark line and a leptonic line, connected by photons. The quark line interacts 
with the background color which is used to describe the high energy nucleus. 
At lowest order in the electro-magnetic and strong coupling constants, only one 
bare photon connects the two fermionic lines, as illustrated in figure (|l|). 

The part of the diagram that describes the scattering of the quark in the 
color field of the nucleus is identical to what we have already calculated for the 
case of photon production. Namely, the photon can be attached before or after 
the scatterings, but the terms where the photon is attached between scatterings 
of the quark on the nucleus are suppressed by inverse powers of the center of 
mass energy yfs. Therefore, compared to the photon production case, we need 
only to replace the photon polarization vectors of the produced photon by the 
propagator of the (now virtual) photon, its coupling to the leptonic line and the 
spinors of the pair. This amounts to the following substitution^ 

*»(feK(fc) - > J+^Jl^ ikuk^) , (4) 



1 We have indicated the ie prescription for the photon propagator, but it is in fact irrelevant 
here since the photon must have an invariant mass squared k 2 larger than Am 2 because of the 
threshold for the production of a pair of leptons with mass rrt; . 
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where k = k\ + It2 is the 4-momentum of the virtual photon, and where L p<T (k) is 
the discontinuity of the one- loop leptonic contribution to the photon polarization 
tensor (i.e. the loop on the upper part of the diagram of figure |l| - only its 
discontinuity is needed since the leptons are produced on-shell). 

An important simplification we used in the photon case was that the sum 
over the photon polarizations turned the product e M (fc)e*(fc) into — g^ v (up to 
terms proportional to k^ that do not contribute thanks to Ward identities). 
This property led to a dramatic simplification of the Dirac algebra involved 
in the calculation of the photon production cross-section. A priori, the object 
L pa {ki : k2) is not proportional to g pa , which means that the fully differential 
lepton pair production cross-section has a more complicated Lorenz structure. If 
however we assume that one reconstructs the virtual photon 4-momentum from 
the momenta k%, k^ of the components of the lepton pair, the same simplification 
occurs. Indeed, if one integrates over the lepton momentum inside the leptonic 
tensor L ptT , keeping the sum k — k\ + fixed, we obtain the following result: 

L p - = ^a cm {g p °k 2 - k p k°) . (5) 

Therefore, for the cross-section da /o I k, we have the same simplification as in 
the photon production case. 

The relation between the differential cross-section and the amplitude can 
therefore be written as 

where M p is the amplitude for the production of a virtual photon, amputated 
of its external legs, and from which the factor 2Tr8(q~ + k~ — p~) has been 
removed. Explicitly, we have: 



M^iplqk) = -ie g u(q) 



1-(tf-$ + m)Y 7 M (</+ # + m)7" 



u(p) 



(p — k) 2 — m 2 (q + k) 2 — m 

d 2 xxe i(g±+k± - p±yx± (U(x±) - l) , (7) 

where e q is the electrical charge of the quark and where U(x±) is a matrix in 
the fundamental representation of SU (N c ) that represents the interactions of 
the quark with the classical color field of the nucleus: 

U(x ± ) = Texp^-ig 2 J dz~ ^-p«(*~ , ^)* a j (8) 

with t a in the fundamental representation, and where p a (z~ , z±) is the density 
of color sources in the nucleus. The color averages over the distribution of hard 
color sources in the nucleus are identical to the case of photon production, and 
can be found in section 3 of 101 . 
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The factor Ad** Ad*, in the cross-section differs from the factor we denoted 
\M\ in our calculation of photon production only in the fact that we must 
not assume k 2 
(tr(£ti)\ 

\ v ' I spin 

value is now: 



in the calculation. In fact, in this quantity, only the factor 
that contains the Dirac algebra is affected by this change. Its new 
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(9) 



where m is the mass of the quark and where we denote D p = (p — k) 2 — m 2 = 
—2p-k + k 2 and D q = (q + k) 2 — m 2 = 2q ■ k + k 2 . It is trivial to check that 
it reduces to the value found in Jj| if we set k 2 = 0. Having in mind the fact 
that the quark comes from the wave-function of a proton, we neglect the mass 
of the quark in the following. 



3 Diffraction 

Like in Q , we can first study the case of diffractive dilepton production, as the 
kinematics is simpler. Let us just remember that the diffractive cross-section is 
obtained by performing the average over the distribution of nuclear color sources 
before squaring the amplitude. This implies that no net transverse momentum 
is exchanged between the nucleus and the quark, i.e. p± = q± + k±. If one 
evaluates the factor (tr(U L)} s ^ with this kinematical constraint, one obtains 
a vanishing result if we neglect the mass m of the quark: 

<^)C - ■ ( 10 ) 

This is in fact similar to the case of real photon production: the absence of 
transverse momentum exchange between the quark and the nucleus prevents 
the emission of the virtual photon. 



4 Inclusive cross-section 

If we do not require a diffractive process on the nuclear side, we just have 
to perform the average over colors sources after squaring the amplitude. The 
details of this procedure are given in section 6 of B . We obtain the following 
expression for the differential cross-section: 

d 4 k d 3 q e 2 q nR 2 2a cm , f . 
dCTincl ~ (2^)4 (2^) 3 2g 2p~ ik 2 L)) ^ 

X2n8{q- +k~ -p~)C(p± q± - k±) , (11) 
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where we define again 



C(I X ) = / d 2 x ± e il ^ (U(0)U*(x ± )) o . (12) 



Like in the case of real photon production, all the information about the nature 
of the medium crossed by the quark (in particular, all the dependence on the 
saturation scale Q s ) is contained in this function C . 

At this point, it is useful to introduce the longitudinal momentum fraction of 
the virtual photon z = kr /p~ , as well as the total transverse momentum transfer 
between the nucleus and the quark l± = q± + k±. The phase space d 4 k of the 
lepton pair can be rewritten as d 4 k — \d(M 2 }{dzj z)d 2 k±, while the phase space 
of the outgoing quark can be written as d 3 q/ (2n) 3 2qo — \(dq~ /q~)8(q~)d 2 l± 
(implicitly, q + = q\j (2q~)). In terms of these new variables, the inclusive 
differential cross-section reads (for unit electric charge of quark) 

it Ft? dz d 2 k± d log M 2 3tt (2tt) 4 1 ±} 



1 +(!-«) 



2;2 



Z 2 l 



-z(l - z)M 2 



[k 2 ± + M 2 (l - z)][(k± - zlx_) 2 + M 2 {1 - z)] 

-l 2 ' 



1 1 



[k 2 ± + M 2 (l - z)] [(k± - zl ± ) 2 + M 2 (l - z)] 



(13) 



This is our main result. It gives the differential cross section for inclusive pro- 
duction of dilcptons in high energy quark- nucleus collisions and includes all the 
high gluon density effects in the nucleus. All the information about the high 
gluon density effects in the nucleus is contained in the function C(l±) || [Io| . 
This function behaves as in the l± ^> Q s (perturbative) region and like 

1/f j_ at l± ~ Q s . For l± <C Q s , it is almost flat. Furthermore, the value of l± 
where the slope of the cross section changes strongly depends on rapidity. This 
slow down happens at higher values of lj_ in the forward rapidity region. 

This expression reduces to the one found in Q for real photon production 
if we take the limit M 2 — ► as it must. The main difference compared to 
the production of a real photon is the fact that the collinear singularities at 
k±_ = zlj_ and k± = are now screened by the invariant mass squared of the 
lepton pair, via the term M 2 (l — z) > 0. 

In order to relate (|l^) to proton-nucleus collisions, we will need to convolute 
( p^[ ) with the quark distribution function in a proton using collinear factorization 
theorem. Explicitly, 

, p A^q l+l-X , . q A^q i+r X 
a °incl ^ / J T „C T r>2N " g incl (-,a\ 

dzd 2 k ± d\ogM 2 J q[ ^ f> dzd 2 k ± d\ogM 2 { ' 

Furthermore, one will need to convolute the above cross section with a quark/ha- 
dron or quark/jet fragmentation function if one is interested in measuring both 
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the outgoing hadron/jet as well as the dilepton. Otherwise, one can do the l± 
integration above to get the pA — ► l + l~ X differential cross section. 

Experimentally, one will be able to study M 2 , k\ and rapidity (z) depen- 
dence of the dilepton production cross section in p — A collisions at RHIC in the 
near future. Here, we outline our qualitative predictions from the Color Glass 
Condensate picture of a nucleus at high energy which will be straightforward to 
verify /falsify at RHIC. 

First, as compared to the standard leading twist perturbative QCD, we 
expect that the partonic level cross section da/dy d 2 l± will change its behavior 
from l/lj_ to for l± ~ Q s and an even flatter behavior at smaller l±, at 
fixed rapidity. Convoluting the partonic cross section with parton structure 
functions in order to get the proton-nucleus cross section will change the power 
of l±. Nevertheless, we expect the difference in the power of l± to be observable 
even after the convolution Q. 

Second, the change of the slope of the cross section from to will 
happen at a higher transverse momentum in the forward rapidity region than 
the mid-rapidity region. Indeed, the saturation scale of the nucleus near the 
fragmentation region of the proton is much larger than the value Q 2 S ~ 1 — 2GeV 2 
usually quoted at mid-rapidity [ pi] ]. The growth of the saturation scale with 
energy is known from DIS experiments at HERA Q and heavy ion collisions at 
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For the same reason as above, transverse momentum broadening of the 
jet+dilepton system will depend on its rapidity: it will be larger at forward 
rapidities. This broadening proportional to Q s adds up to the broadening due 
to initial "intrinsic transverse momentum" of the incoming quark. This effect 
has been neglected in our final result by setting p± = 0. All it would take to 
keep both effects simultaneously would be to keep p± ^ in the calculation. 

A more quantitative investigation of our results is beyond the scope of this 
work and will be pursued elsewhere. Nevertheless, we would like to point out 
that inclusion of the standard leading order pQCD diagrams jt6) such as dilepton 
production via quark-anti-quark annihilation^ and from direct (virtual) photon 
diagrams will be required for numerical accuracy. Since these diagrams are not 
effected by the strong classical field of the nucleus and do not interfere with 
the diagrams considered in this work, one can just add their contribution to 
our results. This way, one will have the full LO [O(of m )] and NLO [0(a s a 2 m ] 
dilepton production cross section in p— A including the high gluon density effects 
in the nucleus 

We would like to emphasize that our result (Jl3|) can also be used for heavy 
ion collisions in the very forward rapidity region where valence quarks are the 
dominant partons in the projectile nucleus. The only difference with p~ A is that 
one would then need to convolute our cross section with the quark distribution 
function in a nucleus rather than a proton. This may make it possible to extract 
the shadowing function for quarks (at large x) in the projectile nucleus by 

2 In principle, one will have to include the effects of high gluon densities on sea quarks. 
However, this will have another factor of a s and therefore is higher order 
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considering the ratio of dilepton cross sections in A — A and p — A collisions at 
RHIC or LHC. 
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